BMC Evolutionary Biology

Early high disparity, where clades reach maximum morphological diversity early in their history, is often observed in the fossil record. Since there are currently no suitable methods to model these Early High Disparity (EHD) processes accurately, we introduce a new phylogenetic comparative method (PCM) to address this gap by combining two evolutionary modes--Ornstein-Uhlenbeck (OU) and Early Burst (EB)--leading to an initial surge in disparity until an early maximum, driven by the EB component, after which disparity declines due to OU constraints. We tested the EHD model on ichthyosaur skull length data and found it to be the best fit in nearly all cases compared to other modes. The median time of peak disparity (TEHD) is about 231.1 million years ago, during the Triassic, roughly 12.4% into the ichthyosaur lineage. To explore how environmental and geological factors might influence peak disparity, we performed changepoint analyses on sea level, sea surface temperature, rock area, and sedimentary sequence ages. While changepoints in temperature at 237.2 Ma and sequence boundaries at 238 Ma were closest to TEHD, no significant changepoints appeared to match TEHD.

The evolution of genome size and eukaryotic genome organisation has puzzled evolutionary biologists for decades. Are differences in DNA content and protein-coding capacity of eukaryotic genomes a result of natural selection, or neutral processes? If they were shaped by natural selection, which features might be connected with these differences: cell volume, ecological strategy, availability of nutrients, etc.? Kinetoplastids are protists with compact nuclear genomes that almost entirely lack introns, but whose mitochondrial DNA content (kDNA) can rival that of their nuclei. We hypothesised that the reduction in non-protein coding nuclear DNA content during the evolution of kinetoplastids is connected with an increase in kDNA content: species with larger fractions of total cellular DNA that is allocated to kDNA will have concomitantly smaller fractions of non-coding DNA in the nuclear genome. We developed a framework for estimation of DNA content in both the mitochondria and nuclei from quantitative fluorescence microscopy data, and inferred the fraction of non-coding DNA in nuclear genomes based on the combination of transcriptomic and fluorescence imaging data. Across 10 genetically distant isolates from diverse free-living kinetoplastid lineages, there is indeed a trend towards a lower proportion of non-protein-coding nuclear DNA when kDNA proportion is higher relative to the total DNA in the cell. This observation is consistent with the hypothesis that in free-living kinetoplastids, mitochondrial and nuclear genomes might co-evolve in terms of their DNA content. This finding raises the possibility that natural selection acts on genome size evolution in kinetoplastids, and contributes evidence and methodology to address this question more broadly in eukaryotes.

The allorecognition system of the hydroid Hydractinia symbiolongicarpus plays an essential role in mediating self vs. non-self recognition among cnidarian colonies, analogous to the role of the major histocompatibility complex (MHC) in bilaterian immune systems. Both complexes have contiguous stretches of similar, highly duplicated genes, suggesting a homologous relationship. Protein structural prediction methods such as AlphaFold2 also indicate structural similarity between the proteins encoded within the Hydractinia ARC and the MHC. We have determined the structure of the ARC in the related hydrozoan Podocoryna americana to determine the degree of conservation of the cnidarian allorecognition complex. Protein structural prediction methods have revealed a number of candidate allorecognition proteins that form a long contiguous complex. Phylogenies suggest homology between cnidarian allorecognition genes, as well as syntenic conservation of this gene complex across hydrozoans. Ongoing comparative genomic analyses across bilaterian species aim to clarify the relationship between the ARC and the MHC, thereby informing our understanding of the evolutionary origins of bilaterian adaptive immunity.

Signals of positive selection around the SLC30A9 gene have been reported in human populations outside Africa. Selection likely acted on a highly differentiated single-nucleotide polymorphism, rs1047626, leading to a non-synonymous substitution in the encoded zinc transporter. Because of the striking similarity between the putatively selected SLC30A9 haplotype observed in several current human populations and the Denisovan individual, previous work has proposed adaptive introgression. Yet alternative explanations, including ancient human variation, and the precise archaic source -Neanderthal or Denisovan- remained unresolved. Considering the potentially complex evolution of SLC30A9, we applied Approximate Bayesian Computation (ABC) algorithms coupled to machine learning to investigate the most plausible evolutionary origin of this substitution. After modelling different evolutionary scenarios with forward-in-time simulations, our results highlight that the most probable scenario is a Denisovan origin of the rs1047626 polymorphism. However, the allele likely introgressed into Neanderthals first and was then passed into non-African modern humans. Moreover, the derived allele frequency for rs1047626 across several African populations is consistent with back-to-Africa migrations. Finally, our ABC analyses indicate strong positive selection in East Asian populations and other out-of-Africa populations, whereas in Neanderthal populations, the selection coefficient was probably neutral or slightly deleterious.

Lineages that have undergone little apparent evolutionary change over long geological intervals are often described as "living fossils", with iconic examples including the coelacanth, tuatara, and Ginkgo. However, studies of morphological and molecular data have produced mixed evidence regarding whether living fossils have evolved more slowly than their relatives. We tested for reduced evolutionary rates using a sister-pairs approach, analysing a comprehensive set of morphological and molecular data for 45 living-fossil lineages and their close relatives across the Tree of Life. Living-fossil lineages showed substantially reduced rates of morphological evolution, whereas evidence for reduced molecular rates was weaker. To test whether this pattern could be explained by species poverty, we also analysed other depauperate lineages (species-poor clades contrasted with species-rich sister clades) and found a weak tendency towards reduced molecular rates but without a corresponding slowdown in morphological evolution. Our study reveals that the living-fossil concept captures a quantifiable pattern of prolonged morphological conservatism, but also suggests that living fossils lack other shared biological or life-history characteristics that might unite them as a group.

The relative roles of genetic drift, inbreeding, and natural selection in island evolution have been fundamental themes in evolutionary biology for nearly a century. However, empirical studies have yet to reach a consensus, and research integrating both genotypes and phenotypes in primates remains limited. The present study aims to estimate the demographic history of an island primate population, the Yakushima macaque (Macaca fuscata yakui), and evaluate the genomic and morphological consequences of its isolation. To this end, we conducted population genetic analyses based on whole-genome sequencing and geometric morphometric analyses of cranial shape. Our results indicate that M. f. yakui has been isolated with an effective population size of several thousand individuals for over 100,000 years. Although we detected signatures of purifying selection, deleterious alleles have not been sufficiently purged. The degree of fluctuating asymmetry in males was slightly, yet significantly, higher in M. f. yakui than in mainland populations, whereas this trend was not observed in females. Marginal evidence of directional selection was found for cranial size, but not for shape. Overall, our findings suggest that due to long-term isolation, genetic drift and inbreeding have exerted a substantial influence on shaping genotypic and phenotypic variation in this island population.

Gene duplication is an important evolutionary mechanism contributing to genomic innovation and adaptation. While the long term fate of duplicated genes has been extensively studied, immediate effects such as those resulting from increase in gene dosage remain understudied. Of particular interest is the question of how expression level relates to the optimum after a duplication, as it could fall short, exceed or directly reach the optimal expression. In this study, we investigated the relationship between expression and fitness for 11 genes previously reported to increase fitness when duplicated in Saccharomyces cerevisiae to compare the expression level achieved by duplication to the one maximizing fitness. While most of the studied genes show little fitness variation across the expression levels tested, we found that ESP1, which encodes a separase essential to anaphase initiation, exhibits a clear fitness optimum at an expression level above the native single-copy state. However, the expression level associated with gene duplication exceeds the optimum. We find that ESP1 overexpression is associated with a reduced cell size and a lower proportion of the cell cycle spent in G2 phase, consistent with its role in cell cycle progression. We propose that the duplication of ESP1 accelerates the cell cycle, resulting in smaller, faster-dividing cells. When measured in different environments, the fitness functions also exhibit similar shapes, with slight variation in the optimal expression. Our findings highlight how gene duplication can confer an immediate selective advantage through dosage increase while being non-optimal, reversing the directionality of natural selection on expression level after duplication.

Random genetic drift, rather than selection, is expected to determine the dynamics of alleles with fitness effects that fall below the reciprocal of effective population size (Ne). This is the drift barrier. The drift barrier hypothesis suggests that adaptation may be hopelessly constrained in organisms with low Ne, with much of the complexity of life reflecting a byproduct of ineffective natural selection. Yet decades of empirical research in evolutionary ecology provides substantial evidence of adaptation in finite metazoan populations. Here I explore the extent to which the drift barrier may, or may not, constrain contemporary adaptive divergence under directional selection in the wild. Compensatory directional selection on polygenic traits occurring after the fixation of a deleterious allele will be twice as efficient compared to the initial action of stabilizing selection. Consistent with this result, meta analysis reveals a negative relationship between contemporary estimates of Ne and the strength of directional selection in natural populations. Analysis of empirical phenotypic selection gradients from wild metazoans further reveals that directional selection is often strong enough to overcome the constraints imposed by finite population size, across a range of genetic architectures. Moreover, on a multi-peak adaptive landscape, the probability of a peak shift between adaptive zones increases exponentially as Ne declines, with observed values of contemporary Ne typically being well within the range needed for peak shift models to be a viable explanation of diversity. Together, these results suggest that widespread adaptation in metazoans is not inconsistent with the fundamental constraints on perfection imposed by drift.

Aquaporin-10 (AQP10) is an aquaglyceroporin that transports small, uncharged molecules, such as water and glycerol. In some species, such as mice and cattle, the AQP10 gene (Aqp10) is pseudogenized. In rodents, pseudogenization of Aqp10 is limited to the Myomorpha suborder, while Aqp10 remains intact in other rodent lineages. The order Cetartiodactyla encompasses diverse species, including Suina, Tylopoda, Ruminantia, and Whippomorpha. The distribution of intact Aqp10 across the entire Cetartiodactyla and analyses estimating the timing of Aqp10 pseudogenization have not been sufficiently investigated. In this study, we identified species within Cetartiodactyla that retain intact Aqp10 and analyzed the detailed patterns of pseudogenization to elucidate the specifics of pseudogenization within this order. We collected Aqp10 loci by BLAST and synteny analyses of publicly available genome databases of 51 cetartiodactyl species from 17 families. Dot plot and alignment analyses revealed that seudogenization of Aqp10 occurred in 42 Cetartiodactyla species, including the Camelidae, warthog (Suidae), Antilocapridae, okapi (Giraffidae), Cervidae, Bovidae, Hippopotamidae, and Cetacea. In contrast, nine species retained an intact Aqp10: pig and babirusa (Suidae), peccary (Tayassuidae), mouse deeres (Tragulidae), giraffes (Giraffidae), and musk deers (Moschidae). The timing of pseudogenization was estimated by comparing patterns of exon deletions, frameshifts, and nonsense mutations in Aqp10 pseudogenes. RT-PCR analysis of cattle and sheep tissues revealed that these species express Aqp3 and Aqp7 in their intestines, but not Aqp10. Collectively, these findings indicate that Aqp10 pseudogenization occurred independently in multiple cetartiodactyl lineages and provide new insights into the evolutionary history of Aqp10 loss within this order.

Many pathogenic and microbial species undergo complex life cycles wherein they experience drastic bottlenecks during invasion into a host and rapid growth within hosts. The effects of such a complex population history on patterns of genetic variation is not currently understood. We model a pathogen population as a metapopulation where each deme represents a host and incorporate transmission dynamics as well as the effects of transmission bottlenecks. We employ a coalescent framework, obtaining analytical expressions for pairwise times to coalescence and genetic differentiation within and between hosts. We find that recurrent bottlenecks rescale the coalescent process within hosts to the size of the bottleneck, reducing within-host variation. In addition, the pairwise distribution of time to coalescence under our model deviates substantially from the Wright-Fisher process when the number of hosts (d) is of similar order to the size of the bottleneck (k), i.e., k [~] d and when m <<1/(2k) +1/(2d), where m is the migration rate. Otherwise, the pairwise process under our model converges in distribution to the Wright-Fisher process. Coalescent simulations show that under the same conditions the site frequency spectrum (SFS) under our model deviates from that under the Kingman coalescent and is drastically skewed towards singletons. Our results suggest that the inference of selection or demographic history from pathogen population genetics data may be systematically biased by the mechanics of pathogen reproduction, indicating a need for further theoretical work that can inform inference.

Mammalian brains exhibit extensive diversity in size, cellular composition, and organization, yet the molecular evolutionary changes associated with this phenotypic diversification remain incompletely understood. Here, we investigated whether protein-sequence evolution covaries with quantitative differences in brain phenotypes across mammals. We integrated 29 brain and body traits with orthologous protein sequences from 18 mammalian species and obtained gene-specific molecular evolutionary measures for 10,395 genes. To reduce redundancy among strongly correlated phenotypes, we selected five representative traits that explained 88.12% of the standardized trait variance under linear reconstruction. Using RERconverge, we tested for associations between gene-specific relative evolutionary rates and evolutionary changes in these five representative traits. The primary analysis, based on log-transformed phenotypes, identified 25 significant gene--trait associations involving 23 genes, 19 of which were also recovered using untransformed phenotypes, indicating that a substantial subset of the associations was robust to phenotype transformation. Expression profiling of these shared candidate genes across adult human GTEx brain tissues revealed heterogeneous patterns ranging from broad expression across the examined tissues to low adult brain expression. Functional enrichment analysis further showed that the gene set retained for comparative analysis was enriched for diverse biological processes, including metabolic, cellular, and neural pathways. These results identify gene-specific evolutionary-rate associations with quantitative mammalian phenotypes and provide functional context for the genes represented in the comparative analysis. Our study provides a comparative framework linking relative evolutionary rates of protein-coding genes to quantitative variation in mammalian brain cellular architecture while accounting for shared evolutionary history.

Recombination shapes allelic combinations across the genome and plays a central role in adaptation and speciation. While recombination rates are known to vary and evolve, how recombination landscapes diverge between populations and which processes drive this divergence remains poorly understood. Here, we use the Population Recombination Divergence Index (PRDI) to quantify divergence in effective recombination landscapes between seven pairs of populations and species from four percomorph fish lineages (cichlids, sticklebacks, killifish, and clownfish). While accounting for stochasticity and confounders of linkage disequilibrium (LD) -based recombination rate estimates, we find that recombination landscapes diverged substantially in all comparisons. This divergence generally increases with genetic differentiation and divergence time, but the amount of PRDI accumulated per unit of genetic divergence (FST), net nucleotide divergence (da), or divergence time varies markedly among population and species pairs. Comparisons across systems suggest that PRDI accumulation tends to be faster (per unit of FST, da, or divergence time) in pairs with higher inferred gene flow. We also introduce a local version of PRDI (_l) to better characterize genomic heterogeneity in effective recombination divergence. We find that higher inferred gene flow between populations tends to be associated with stronger intragenomic correlation between local FST and _l. Finally, we propose a conceptual framework for investigating how divergent selection and gene flow may interact to shape effective recombination landscape divergence. Our work highlights the heterogeneous dynamics of effective recombination landscape evolution, and motivates further theoretical and empirical efforts to integrate recombination landscape evolution into studies of adaptation and speciation.

Fossil central nervous systems (CNS) have illuminated several facets of total-group euarthropod evolution, including their phylogenetic affinities, the origin of the head, and their ancestral neuroanatomical organization. However, the CNS remains unknown in trilobites, the most diverse and earliest group of Palaeozoic euarthropods, but can potentially clarify fundamental questions that remain disputed despite thousands of described species and countless studies spanning over two centuries, such as their affinities relative to modern representatives. Here, we describe the CNS of the cheirurid Anacheirurus adserai from the Lower Ordovician Fezouata Shale. Anacheirurus features a brain with well-developed proto, deuto- and tritocerebral neuromeres, an anteriomedian protocerebral protrusion, a central complex, an esophageal foramen, and a ganglionated ventral nerve cord. The CNS is replicated in framboidal iron that differs in texture and composition from the exoskeleton and rock matrix. The overall organization of Anacheirurus's CNS most closely resembles that of Cambrian stem-group mandibulates, and critically shares the presence of a central complex with a protocerebral bridge and lateral accessory lobes with extant pancrustaceans. Our findings provide new morphological support for the hypothesis that Trilobita belongs to the mandibulate stem lineage and bring new insights into the ancestral morphology and evolutionary history of the euarthropod central complex. More broadly, our findings offer direct anatomical evidence of the key sensorimotor system innovations that originated during the Cambrian Information Revolution.

Ancestral recombination graphs (ARGs) represent the coalescent history at every locus in the genome, subject to mutation and recombination. Modelling viruses using the coalescent with recombination remains less explored, due to the potential for complex, multi-scale dynamics. Using perturbative methods, we establish mathematical bounds that describe the balance between inter-host and within-host coalescence, and demonstrate that the highly-recombinant Epstein-Barr Virus (EBV) coalesces predominantly at the inter-host level. We then modify SINGER, a tool from diploid genetics, for haploid ARG inference with fine-grained variable recombination rates. We provide insight into the paradox of high recombination rates and slow linkage disequilibrium decay, as ARG-based methods potentially allow us to disentangle epistatic effects from shared ancestry, and illustrate this using the coalescent topologies of EBV's host cell-entry glycoproteins. Our findings demonstrate that ARG-based methods recover gene-level selective structure, and provide a general method for analysing virus phylodynamics in the limit of small quantities.

Whether sexual selection facilitates or hinders adaptation to novel environments is a question of fundamental importance due to current rapid environmental change. Using phylogenetic comparative analysis in 42 Drosophila species, we found that sperm length - a sexually selected ornament - positively associates with the ability of males to survive high temperatures. However, the ability to remain fertile does not increase with sperm length to the same extent, meaning that longer sperm species experience a larger gap between temperatures that kill and temperatures that sterilize males. The findings are consistent with the idea that strong sexual selection on condition-dependent exaggerated traits may facilitate somatic resilience to environmental stressors, but that susceptibility to heat-induced infertility may be a cost of extreme sperm ornamentation. Overall, the data suggests multifaceted links between postcopulatory sexual selection and climate change.

Single-cell sequencing provides a rich source of data that can be used to model and estimate the evolution of somatic cells within individuals. Many common models for modeling DNA evolution in phylogenetics are based on continuous-time Markov processes that operate on aligned sites. For haploid DNA, such as found in mitochondria, there are many well established models for these processes but since nuclear DNA is diploid, containing copies of both the maternal and paternal genome, it is natural to develop continuous-time Markov processes adapted to diploid DNA. One possible model is that both the maternal and paternal evolve independently according to haploid models so that, in principle, DNA comes with a phase. i.e. a label that tells if the site was derived from the maternal or paternal sequence. However, this information may not be available in practice and only an unphased diploid genotype is available. That is, an unphased model groups both of the possible maternal-paternal labels into a single state. This grouping process, called ``lumping'', is used by other methods as well and results in a coarser process. Care must be taken so that the coarser process can still be modeled as a continuous-time Markov chain; the lumped process of the original phased Markov chain may or may not be Markovian, depending on restrictions imposed on the generator matrix. Here, we explore the mathematical formulations of phased and unphased Markov processes, as well as a collection of other diploid DNA models from the literature.

Despite sharing nearly identical enzymes, strains of the same microbial species display wide variability of growth rates when growing on the same nutrients. Here, we show that such growth rate variability can emerge as a consequence of evolutionary adaptation to different patterns of environmental fluctuations. We develop a mathematical model combining cellular proteome allocation with eco-evolutionary dynamics in fluctuating environments. We find that different patterns of environmental fluctuations select for distinct growth profiles. Slow-growing strains repeatedly outcompete fast-growers when evolved in rapidly changing environments. In a range of environments with asymmetric nutrient availability, evolution reproducibly leads to the robust diversification of two coexisting strains. We develop a theoretical framework based on adaptive dynamics which quantitatively predicts these evolutionary outcomes and explains them in terms of the growth-lag tradeoff that all cells face. Finally, we show that the model can reproduce several patterns in growth-rate data from closely-related strains. Our results highlight that bacterial growth rates can rapidly adapt to their recent environmental history by proteome allocation alone, in ways that often oppose expectations based on nutrient quality.

Demographic inference based on sequentially Markovian coalescent (SMC) model is widely used for disentangling the evolutionary history of species and populations. However, the genetic effects of self-fertilization, including the expected reduction of population-level genetic diversity and hence effective population size (Ne), are often ignored in demographic inference, despite the fact that selfing is a common reproductive mode in both invertebrates and plants. To investigate how self-fertilization could bias demographic inference, we apply eSMC2, a selfing-aware demographic inference method, to Hawaiian populations of Caenorhabditis elegans, a predominantly self-fertilizing nematode. Using an expanded dataset that contains more Hawaiian strains than in previous studies, we recovered previously reported genetic groups, including one largely consisting of strains from the Big Island. Population substructure is also detected among these Big Island subpopulations, highlighting the importance of identifying subgroups to avoid misinterpretations of demography. For different Big Island subpopulations, we estimate high selfing rates compatible with previous studies and infer consistent population declines. Population sizes inferred using eSMC2 are smaller than those inferred using PSMC', which does not account for selfing. However, rescaling PSMC' parameters to account for self-fertilization produces a similar demographic trajectory as eSMC2 inference. Our study thus demonstrates how one can account for effects of selfing when using SMC-based methods. It also indicates that the Hawaiian group shows population declines similar to those previously found outside Hawaii.

Crozier's Paradox states that genetic kin recognition will select for its own demise, and theoretical analyses over the past few decades have supported this conclusion. However, empirical evidence is far more limited. Here we examine support for Crozier's Paradox in two kin recognition genes in the social amoeba Dictyostelium discoideum, tgrB1 and tgrC1, which enable cells to recognize and reject unrelated cells during cooperative multicellular development. We find that co-isolated amoebae, taken from just a few grams of soil, harbor highly diverged recognition alleles, indicating plentiful sequence variation for social decision-making on a micro-scale. tgrB1 and tgrC1 belong to a large gene family, whose members vary in their levels of polymorphism, strength of selection, and retention of polymorphism. We do not see strong evidence for the paradox, in the form of positive frequency-dependent selection, or its outcome, in the form of depleted genetic variation. Consistent with previous findings, we detect recombination among the duplicated tgrB and tgrC genes, which suggests a mechanism to generate novel alleles and replenish genetic variation. We suggest that the evolutionary dynamics of large gene families may be essential for rapidly generating novel variation that can prevent Crozier's Paradox from taking effect.

Theory predicts that populations should evolve mutational robustness the ability to maintain a constant phenotype despite mutations. However, selection does not operate on individual mutations but on phenotypes resulting from entire replication events, which may involve a variable number of mutations. We therefore distinguish mutational robustness from replicative robustness, defined as the probability that a replication event produces similarly fit offspring. Mutational and replicative robustness may seem closely aligned at first sight, but we show with an evolutionary simulation model that they can be selected in opposite directions. Our results show that while genomes with high redundancy and/or large amounts of non-coding DNA are mutationally robust, they are replicatively fragile because their size increases the number of mutations per replication. This is especially relevant when genomes are at risk of severe disruption from chromosomal rearrangements. Our results expose limits of the fitness landscape metaphor in the context of robustness evolution: neutral plateaus may be mutationally robust but replicatively unstable and therefore disfavoured by selection, whereas narrow peaks arising from densely encoded genomes can replicate more reliably. Our results identify replicative robustness as a distinct target of selection that depends not only on tolerance to individual mutations but also on their overall likelihood of occurring.